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ABSTRACT 

Hard xray data from the RXTE observatory (HEXTE energy range 15 to 240 keV) have 
been analyzed to obtain a phase coherent timing solution for the Crab pulsar glitch of 15 July 
2000. The results are: (1) step change in the rotation frequency i/q of the Crab pulsar at the 
epoch of the glitch is Ai/q = (30 ± 3) x 10“® x i/q, (2) step change in its time derivative is Az/q = 
(4.8±0.6) X 10“^ X z>o, and (3) the time scale of decay of the the step change is Td = 4.7±0.5 days. 
The first two results are consistent with those obtained at radio frequencies by the Jodrell Bank 
observatory. The last result has not been quoted in the literature, but could be an underestimate 
due to lack of observations very close to the glitch epoch. By comparing with the monthly timing 
ephemeris published by the Jodrell group for the Crab pulsar, the time delay between the main 
peaks of the hard xray and radio pulse profiles is estimated to be +411 ± 167 fisec. Although 
this number is not very significant, it is consistent with the number derived for the 2 to 16 keV 
energy range, using the PCA instrument of RXTE. The separation between the two peaks of the 
integrated pulse profile of the Crab pulsar, and the ratio of their intensities, both are statistically 
similar before and after the glitch. The dead time corrected integrated photon flux within the 
integrated pulse profile appears to decrease after the glitch, although this is not a statistically 
strong result. This work achieves what can be considered to be almost absolute timing analysis 
of the Crab pulsar hard xray data. 


Subject headings: pulsars: individual (Crab Pulsar) X-rays: stars 


1. Introduction 

ditches in pulsars are events in which the 
rotation period and its derivative (and possibly 
higher derivatives as well) undergo an abrupt 
change in value, on time scales of less than 
minutes, often followed by a recovery to ap¬ 
proximately the pre glitch values, over time 
scales of day s to tens of days, or even much 
longer: see (IShemar fc Lvne _ 1R96 Lvne et al 


2006 : Wong et al.l 2001 : Espinoza et al. 201l[l for 
details of pulsar glitches, their history and their 
relevance. Glitches are important to study be¬ 
cause they are probably one of the very few 
methods available to study the intern al struc¬ 
ture of neutron stars ( Bavm et al.lll969tl : see also 
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( Ruderman et al. 1998[l and references therein. 
Glitches are rare events; in the Crab pulsar (PSR 
B0531+21 or J0534+2200) they o ccur once in 


1.6 years (lEspinoza et al.l 120111) . A typical 


glitch in the Crab pulsar involves a very small 
fractional change of rotation period (or alter- 
nately, rotation frequency) of about 10 ~'^ to 10“® 
( Wong et al.ll200ll : lEspinoza et al. 2011 1. Coupled 
with the sudden onset of a glitch, this implies 
that frequent and regular pulsar timing observa¬ 
tions are required to study the glitch phenomenon. 
The very low rotation period of the Crab pulsar 
{P « 33.5 ms, or alternately very high rotation 
frequency i' ~ 29.851 Hz, at the middle of the year 
2000) necessitates timing observations at least 
twice a day, over a pe riod of one year, to properly 
analyze a glitch; s ee ( Manches ter fc Taylor 1977t 


iBacker fc Helling^ 19861 : Lvne fc SmithI 2006h for 

pedagogical reviews of pulsar timing in general, 
and analysis of pulsar glitches in particular. 
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Clearly a dedicated telescope is required to 
study pulsar glitches. At radio frequencies this 
has been done (and continues to be done) for 
the Crab pulsar b y the Jodrell Bank observatory 
( Lvne et al. 1993h . Over the last ~ 40 years, 
they have accumulated pulse timing information 
of the Crab pulsar at 610 and 1400 MHz ra¬ 
dio frequencies, observing it daily, and have pub¬ 
lished (and continue to update) the m onthly tim¬ 
ing ep hemeri jil of the Cra b pulsar: see (iLvne et ah 
1993li for details; also see ( Lvne et al. Il2015ll for the 


45 year rotation history of the Crab pulsar. The 
Crab pulsar has also been observed daily by the 
Green Bank Telescope, at 327 and 610 MHz radi o 
frequencies ( Backer et al ]l2000l:IWong et al.ll200lh . 
Ideally such work should have been carried out at 
xray energies, which are not affected by problems 
associated with propagation through the interstel- 
lar medium, that radio signals are su sceptible to 
( Lvne et al. 1993L iBacker et al.ll2000ll . However, 
xray telescopes are difficult to build and expen¬ 
sive, in comparison to radio telescopes. RXTE is 
one of the few xray observatories that can time the 
arrival of pulses from pulsars to the accuracy re¬ 
quired for timing analysis. However RXTE is not 
a dedicated pulsar timing observatory, and most 
Crab pulsar observations of RXTE are spaced, on 
the average, two weeks apart. Eortunately, dur¬ 
ing the period late 1999 to late 2000, three sets 
of very closely spaced RXTE observations of the 
Crab pulsar were available, one of them being just 
after glitch of 15 July 2000. It is mainly these 
three clusters of observations, and the existence 
of observations immediately after the glitch, that 
have motivated this work. 

This work is organized as follows. Section [5] and 
section|3]describe HEXTE data and the method of 
analysis. Section |4] presents phase coherent timing 
results for the July 2000 glitch of the Crab pulsar. 
Section [5] contains discussion and describes the be¬ 
havior of some Crab pulsar parameters before and 
after the glitch. 


2. Observations 


The HEXTE instrument ( Rothschild et al 


199811 of RXTE consists of two independent clus¬ 


ters of detectors, labeled clusters 0 and 1. Each 
cluster contains four NaI(Tl)/CsI(Na) phoswich 


^http: //www. jb.man.ac.uk/pulsar/crab.html 


scintillation photon counters, and has a field of 
view of one degree in the sky. Eor practical pur¬ 
poses this instrument is sensitive to photons in 
the 15 to 240 keV range, and each photon’s ar¬ 
rival time is measured with an accuracy of ~ 7.6 
p,sec (see “The ABC of XTE” guide on the RXTE 
websit^. During normal operation, the two clus¬ 
ters switch between the source and a background 
region of the sky, such that when one cluster is 
pointed at the source, the other is pointed at the 
background, and vice versa. For pulsar observa¬ 
tions another mode of observation is also used, in 
which both clusters dwell only on the source. The 
data used in this work consist of both modes, but 
predominantly of the latter kind. 


The first observation used in this work was 
obtained on 18 Dec 1999, and the last on 24 
Dec 2000; the corresponding observation identifi¬ 
cation numbers (ObsID) for the data are 40090- 
01-01-00 and 50804-01-14-00, respectively. The 
epoch of the glitch is MJD 51740.656 ± 0.002 


(lEspinoza et al.l 120111) . which is at Ri 15:45 UTC 


on 15 July 2000. The pre glitch data is relatively 
more frequently observed during the months Dec 
1999 and Jan 2000, but not later. It extends 
up to 14 May 2000 only (ObsID 50099-01-26-00), 
which implies that no observations exist for the 
two month duration just prior to the glitch. For¬ 
tuitously, the first post glitch observation is on 
17 July 2000 at Ri 00:45 UTC (ObsID 50098-01- 
01-00), which is just R:: 1.4 days after the glitch. 
From then onward the data is well sampled (fre¬ 
quently observed) up to 31 Jul 2000 (ObsID 50099- 
01-02-00), after which the data is under sampled 
until 5 Dec 2000 (ObsID 50099-01-11-00), which is 
for most of the post glitch duration. Then again 
the data is well sampled until 24 Dec 2000 (Ob¬ 
sID 50804-01-14-00). There are 50 ObsID during 
this period, out of which 2 were not useful; during 
ObsID 50100-01-01-05 the Crab pulsar was com¬ 
pletely occulted by the Earth (ELV < 0°, which is 
the instantaneous angle between the Earth’s limb 
and the astronomical source); and ObsID 50099- 
01-05-00F was obtained on 11 Sept 2000, during 
the week when there were up to 1° errors in the 
spacecraft attitude. Six of the remaining 48 Ob- 
sIDs have more than 10% data gaps in them, but 
they have sufficient useful data for our purpose. 


heasarc.gsfc.nasa.gov/docs/xte/data_ analysis.html 
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Thus the Crab pulsar is under sampled (from the 
timing point of view) for most of the one year du¬ 
ration under consideration, except for three brief 
periods of well sampled data, one at the very be¬ 
ginning of the pre glitch duration, one just after 
the glitch, and one at the very end of the observa¬ 
tions used in this work. 


In fact the pre and post glitch duration were 
chosen based partly on the availability of closely 
spaced (frequently observed) data. The other rea¬ 
son was for the data to be sufficiently isolated in 
time from the previous and the next Crab pul¬ 
sar glitches, so that the timing analysis is not cor¬ 
rupted by their residual effects. The previous Crab 
pulsar glitch (a small glitch) occurred at « 00:29 
UTC on 1 Oct 1999, and the next glitch (a large 
one) occurred a t « 01:44 UTC on 24 Jun 2001 
( Espinoza et al. 201 ih : the very small glitch at « 
18:00 UTC on 17 Sept 2000 is ignored because 
the data of this work is not sensitive to it. The 
first observation of this work is 78 days after the 
previous glitch, which is several times the longest 
of the short decay timescales for the Crab pulsar 
(iLvne et al.l l200Clt Wong et al. 200 ll: IWang et~aI1 
20121) . so one expects that the pre glitch timing 


solution in this work is not corrupted by the de¬ 
cay phenomenon of the previous glitch. The last 
observation of this work is six months before the 
next big Crab pulsar glitch. 


3. Data Processing 

All data used in this work have been ac¬ 
quired in the Event List mode (operating 
modes E_8us_256_DX0F, E_8us_256T)XlF, etc.), 
in which the photon arrival times have accuracy 
~ 7.6 iisec (the best possible for HEXTE), and 
also the highest energy resolution (256 channels 
in the energy range 0 to 250 keV); see “Reduc¬ 
tion and Analysis of HEXTE data” on the RXTE 
websit^. Throughout this analysis, data of each 
cluster are analyzed separately. 

The first step in data processing is the creation 
of the so called Good Time Intervals (GTl), which 
are time duration identifying the useful data, not 
corrupted by instrumental and extraneous factors. 
The GTl that account for system provided checks 
are created by the tool maketime using the fil¬ 


® heasarc.gsfc.nasa.gov/docs/xte/recipes/cook_book.html 


ter file available in the directory STDPROD, and 
using the screening criterion (1) ELV > 10°, (2) 
the difference between the source position and the 
pointing of the satellite (OFFSET) < 0.02°, (3) 
the time since the peak of the last South Atlantic 
Anomaly passage (T1ME_S1NCE_SAA) > 30 min 
or T1ME_S1NCE_SAA < 0 min. 

3.1. Standard Processing of HEXTE Data 

The next step is to create additional GTl based 
on the light curves of the observation; for exam¬ 
ple, abrupt changes of large magnitude in photon 
count rates should be excluded from the analysis. 

For that, each data file for each cluster is pro¬ 
cessed by the tool hxtback to separate the data 
pertaining to source and sky background regions; 
for timing analysis only the source data are used. 
Then light curves are obtained using the tool se- 
extrct, screening the data using the GTl avail¬ 
able in each file, as well as the GTl file created in 
the previous section. The light curves are binned 
at the telemetry interval DELTAT (16 sec), and 
photons are selected from energy channels 15 to 
240 keV. The light curves are corrected for dead 
time using the tool hxtdead, using the appropri¬ 
ate house keeping file. Although all four detec¬ 
tors are chosen for both clusters while running the 
tools, the third detector of cluster 1 lost ability 
to assign energy information to a photon after 6 
March 1996; photons of this detector fall in the 
first two energy channels irrespective of their ac¬ 
tual energy. Therefore our choice of 15 to 240 keV 
energy range essentially filters out photons from 
this detector, even though they have valid arrival 
time information; see “The XTE Technical Ap¬ 
pendix” on RXTE websit43- 

Light curves for both clusters are plotted, and 
the range of count rates, within which the data 
appears good, are chosen. These light curves and 
count rate limits are used to create a second set 
of GTl files, one for each cluster, using the tool 
maketime. These GTl are then merged with 
the earlier GTl, using the tool mgtime with the 
AND option, to yield the final GTl, which are then 
used to filter the individual data files, using the 
tool fselect. By this stage, one has screened the 
data for all system provided checks, as well as for 
user provided count rate limits. Finally, the pho- 


heasarc.gsfc.nasa.gov/docs/xte/appendix_f.html 
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ton arrival times are referred to the solar system 
barycenter using the tool fexbary, using the or¬ 
bit file for the given ObsID, and the Crab pulsar’s 
coordinates (83.6332208° for right ascension and 
22.0144611° for declination, for the epoch J2000, 
( McNamara 1 ll97lll . 1,aken from the online pulsar 
catalog of ATNeH) . 

3.2. Obtaining the Period of Crab Pulsar 

The next step is to obtain the best period of 
the Crab pulsar for each data file of each clus¬ 
ter. The fundamental frequency in the power 
spectrum of the data gives the first approxima¬ 
tion to the period, which is obtained using the 
tool powspec, with bin size 0.67 ms and data 
length 2^° = 1048576 bins for most files, but half 
or quarter of that for shorter files. This is done 
for each data hie of each cluster, not only to check 
for consistency of the period among all hies of a 
single ObsID, but also to check the health of the 
data. The second approximation to the period is 
obtained by folding the data over a range of 600 
periods centered on the hrst approximation pe¬ 
riod, and searching for the maximum using 
the tool efsearch, with an increment of 10“® sec 
in period, and zero period derivative. The third 
approximation to the period is obtained by do¬ 
ing a hner search centered on the second approx¬ 
imation period, using an increment of 0.2 x 10“® 
sec in period, and a nominal period derivative of 
420 X 10“^® sec per sec (explained later), over a 
range of 128 periods. A Gaussian is ht to the as 
a function of period, to obtain the centroid. The 
hnal period is obtained by folding the earlier and 
later portions of the data at the third approxima¬ 
tion period, then cross correlating the two inte¬ 
grated profiles, then measuring the shift (if at all) 
versus time between the two profiles. Folding is 
done using the tool efold, while cross correlation 
is done using independently developed software. 
The accuracy of the final period was typically 7 
nano sec; for ObsID with long duration observa¬ 
tions it could be as small as fraction of a nano 
sec. 


5 WWW. atnf.csiro.au / research / pulsar / psrcat / 


3.3. Obtaining the Epoch 
peak of Crab Pulsar 


of the Main 


The next step is to combine data files of each 
ObsID to get the integrated profile of the Crab 
pulsar separately for each cluster, with a resolu¬ 
tion of 128 bins per period, using the best period 
of the epoch, and a nominal period derivative of 
420 X 10“^® sec per sec, that is obtained by fit¬ 
ting the periods versus epochs for the 48 ObsIDs; 
this value is also consistent with the mean period 
derivative of the Crab pulsar for the duration Dec 
1999 to Dec 2000, as estimated from radio data, 
given in the Jodrell monthly ephemeris. The zero 
phase of the integrated profile is set to the ar¬ 
rival time of the first valid photon in the data, af¬ 
ter conversion to MJD by adding the MJDREF 
available in the data files. For validity of this 
procedure, it must be ensured that the keyword 
RADECSYS in the barycenter corrected data files 
is set to EK5, the keyword TIMEZERO is set to 
0, and the keyword CLOCKAPP is set to T; see 
“A Time Tutorial” in “The ABC of XTE” guide, 
and also “RXTE Absolute Timing Accuracy” on 
RXTE website0. 

The final step is to obtain the epoch of arrival of 
the peak of the main pulse of Crab pulsar, which is 
taken as the fiducial point in its integrated profile. 
For this the epoch of the zero phase of the inte¬ 
grated profile (described above) should be added 
to the position of the peak of the main pulse. This 
position is found by three independent means, as 


done in ( Rots et ^ 2006ll but with some differ¬ 


ence - (1) fitting a Gaussian to the main peak 
data, (2), fitting a Lorentzian to the main peak 
data, and (3) finding the first moment of the main 
peak data higher th an 80% of the pea k value. The 
first difference with (IRots et al.ll2006l) is that they 
fit a parabola instead of a Gaussian in method (1). 
The second difference is that they have 200, 400 
and 800 bins in the integrated profile for the three 
methods, respectively, whereas the number of bins 
in this work are 128 for all three methods, because 
several ObsID do not have sufficient exposure to 
obtain sufficient signal to noise ratio with higher 
number of bins. The disadvantage of having lower 
number of bins in the integrated profile is that 
one has to be cautious while fitting the Gaussian 
and the Lorentzian; one has to choose as much of 


® heasarc.gsfc.nasa.gov/docs/xte/abc/time.html 
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the main peak data as possible, to maximize the 
sensitivity of the fit, but should not include the 
asymmetric parts of the peak in its wings. This 
fitting is done for each cluster for each ObsID. All 
three methods give consistent timing results for 
the Crab pulsar glitch of July 2000. The typical 
accuracy of the fit is Ri 1 milli period. Figure [T] 
shows the integrated profile of the Crab pulsar for 
ObsID 40090-01-01-00 for the combined data of 
both clusters. 



Fig. 1.— Integrated profile of Crab pulsar obtained 
by combining data of ObsID 40090-01-01-00 from both 
clusters. The epoch of the peak of the main pulse of 
cluster 0 data is used as reference phase for the tool 
efold, and data is folded at the period of this epoch 
(0.0335046788 sec). The average counts per sec (used 
to normalize the ordinate) is 176.477783. 

4. Phase Coherent Timing Solution for 
the July 2000 Glitch 

The timing solution for Crab pulsar presented 
in this work was obtained using the Gaussian fit¬ 
ting method described earlier, with three vari¬ 
ations of epoch of arrival of the peak of the 
main pulse (henceforth referred to as pulse ar¬ 
rival epoch) - (1) epochs from data of cluster 0 
only, (2) epochs from data of both clusters, but 
separately, and (3) epochs from combined data of 
both clusters. Methods (1) and (3) both have 48 
epochs, while method (2) has twice the number. 
The epochs from method (3) are expected to be 
the most reliable due to enhanced number of pho¬ 
tons in the integrated profile by the average factor 
1.75 (cluster 1 gathers 25% less photons due to ex¬ 
cluding detector 3). All three methods give consis¬ 


tently similar results. The results presented in this 
work are derived using method (3). As mentioned 
earlier, timing solutions obtained using method 
(3) along with fitting a Lorentzian, and finding the 
first moment, all give results similar to those pre¬ 
sented in Table [2] in this work. In addition, com¬ 
pared to the Gaussian fit method, the first moment 
method gave a mean departure of 1 ± 2 milli peri¬ 
ods for the 48 epochs; for the Lorentzian method 
the mean departure was 0.3 ± 0.4 milli periods. 
Thus it is concluded that pulse arrival epochs are 
consistent among the three peak finding methods. 
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Fig. 2.— Result of using TEMP02 on the 14 pre 
glitch pulse arrival epochs, fitting for v, v and v, 
using the epoch of the first data at -210.221 (MJD 
51530.4349129857) as the reference epoch for phase 
zero in TEMP02; the results are given in Table [T] 
The origin of the abscissa is at the epoch of the glitch 
(MJD 51740.656). Note the four closely spaced epochs 
at the start of the data, and the lack of any observation 
for two months prior to the glitch. 


Table 1: TEMP02 best fit parameters to the pre 
glitch data of Figured] v is the rotation frequency of 
the Crab pulsar at the epoch MJD 51530.4349129857, 
which is the first epoch in our data; v and v are the 
first and second time derivatives of v, respectively, at 
the same epoch. The errors in brackets are in the last 
digit of each result. 


Parameter 

Value 

V (Hz) 

29.846592902(2) 

V (10-1*^ s-^) 

-3.745962(8) 

V (10-^*^ s-^) 

0.94(1) 


Figure [2] shows the result of using TEMP 02 
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( Hobbs et al. 2006f) on the 14 pre glitch pulse ar¬ 
rival epochs. In the absence of glitches and tim¬ 
ing noise, pulsars obey a simple slowdown model 
that is adequately represented, at least for short 
duration of several months, by three parame¬ 
ters - rotation frequency v, its time derivative 
v and its second deri vative v\ see Equation 1 of 
( Espinoza et al ] l2niih . The three fitted parame¬ 
ters are shown in Table [l] and are consistent with 
the interpolated values from the Jodrell monthly 
ephemeris for the Crab pulsar for that epoch. The 
formal one standard deviation error on the pulse 
arrival epoch in Figure [2] is typically 0.87 milli 
periods. An independent estimate of the error is 
obtained by the difference in pulse arrival epochs 
for clusters 0 and 1 for the same ObsID. Often 
these differ by « 16 sec or a few multiples of it, 
but on rare occasions can differ by hours. Ideally 
these relatively short duration differences should 
be equivalent to almost integer number of pulse 
cycles, since the estimated period of rotation of 
the Crab pulsar would be very accurate for closely 
spaced epochs; the departure from integer values 
will give us an idea of the errors involved in pulse 
arrival epochs in Figure [2j The mean value of the 
departure from integer number of cycles turns out 
to be 1.0 milli period. Therefore all error bars 
in Figure [2] are set to one milli period. The rms 
residual of the 14 pre glitch arrival times after the 
TEMP02 fit is 3.2 milli periods. 

TEMP02 is used with the parameters of Ta¬ 
ble [T] as constant input (i.e., without any fitting) 
for the 34 post glitch pulse arrival epochs; the re¬ 
sults are shown in Figure [H These post glitch 
residuals A;/) (in sec) are fit to a modified ve rsion 
of the glitch model of ( Shemar fc Lvnd 1996ll . 


/\4) 



Aup + Avpt + Ai/„ exp 


-A(1)q - 
TdAVri 

va 


Avp ^ Ai>p 

Uq Uq 2 




( 1 ) 


where t is the time elapsed since the glitch epoch 
(in seconds), vq is the rotation frequency of the 
Crab pulsar at the epoch of the glitch (in Hz), 
Avp and Avn are the permanent and exponen¬ 
tially decaying parts of the step change in rota¬ 


tion frequency at the epoch of the glitch, Az>p is 
the permanent step change in the time derivative 
of the rotation frequency, and is the decay time 
scale (in sec) of the step frequency change. The 
parameter At/fp accounts for any uncertainty e in 
the epoch of the glitch, which is assumed to be 
much smaller than Td- The equation abov e differs 
from Equation 1 of ( Shemar fc Lvnd 1996h in hav¬ 
ing a single exponential only; as mentioned in their 
paper, only occasionally one requires more than 
one transient component, and the decay times of 
these additional transients are typically hundreds 
of days. The negative sign in Equation 1 accounts 
for the fact that the phase residuals after a glitch 
increase in the negative directio n for a positive 
step c hange in rotation frequency ( Shemar fc Lvnd 


1996ll . Table [2] gives the minimum fit values of 


the parameters in Equation 1. The rms residual 
of the 34 post glitch arrival times after the fit is 
8.6 milli periods. 
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Fig. 3.— Result of using TEMP02 on the 34 post 
glitch pulse arrival epochs, using the values of Table [T] 
as input parameters and without fitting. The origin 
of the abscissa is at the epoch of the glitch (MJD 
51740.656). The vertical dashed line is at epoch after 
which one phase cycle was subtracted for ah subse¬ 
quent epochs, using the command PHASE —1 in the 
input file to TEMP02. Note the clusters of closely 
spaced epochs, one just after the glitch epoch, and the 
other at the end of the data. 


The epoch of the glitch is set to the published 
value of MJD 51740.65 6, since our data is un able 
to verify this number. ( Espinoza et al. HmH) ob¬ 
tain this by using trial glitch epochs as reference 
epochs in TEMP02 for the pre and post glitch 
timing solutions, then comparing the two solutions 
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Table 2: Minimum parameters obtained by fitting 
Equation 1 to the post glitch data of Figure [21 


Parameter 

Value 

A(/o (ms) 

3.3 ±0.5 

Avp (10-° Hz) 

0.180 ±0.003 

Avp (10-33 s-"") 

-0.350 ±0.006 

Aun (10-° Hz) 

0.71 ±0.08 

Td (days) 

4.7 ±0.5 


for a match in phase at the reference epoch (glitch 
epoch). This is not possible in this work because 
of under sampling of data; the one and only pre 
glitch solution available (TablejT]) was with the ref¬ 
erence epoch used, because close to that the data 
was well sampled. An alternate method of formu¬ 
lating the parameter is to take the integral in 
Equation 1 from the limits 0 to t — e, instead of 
from e to t. It can be shown that the two methods 
are equivalent as long as e/r^ ^ 1, and the derived 
parameters Avp, Avp and have the relative 
orders of magnitude as derived in Table 2. 


The total step change in rotation frequency at 
the glitch, as a fraction of the pre glitch frequency, 
is (Ai/p -|- Ai/„)/i/o = (30 ± 3) X 10“®, which is not 
too different fr om the value of (25.1 ± 0.3) x 10“® 
published by ( Espinoza et al.lEoil b The actual 
step change of Avp + Ai/„ rs 0.89 ± 0.08 /rHz 
compares well with t he top panel of Figure 5 of 
( Espinoza et al.l[201lll . The fraction of frequency 
recovery Q is 0.71/(0.71 -b 0.18) Ri 0.80 ± 0.11, 
which is very high, a s evident from Figure 5 of 
( Espinoza et al ] l201lll , although the decay of fre¬ 
quency in their figure (top panel) does not appear 
to be entirely exponential, most probably due to 
the very small glitch that has been ignored in this 
work. 


The total step in frequency derivative at the 
glitch as a fraction of the pre glitch frequency 
derivative is {Avp — Ai/„/(t * 86400))/z>o = (4.8 ± 
0.6) X10“^, which is also no t too different from th e 
value of (2.9±0.1) x 10“^ of ( Espinoza et akllioilh . 
The fraction of recovery of the frequency deriva¬ 
tive is -17.48/(-17.48-0.35) « 0. 98T0.16, which 
is co nsistent with Figure 5 of (lEspinoza et al. 
I 2 OIIII . However, the step change in i>o « (—18 ± 
2) X 10“^^ estimated here is in consistent wi th 


the b ottom panel of Figure 5 of ([Espinoza et al 
I 2 OIII I. in which it is more like Rs —5 x 10“^^ 


However, the value expected from their work is 
Rs 2.9 X 10-3 X -3.744255 x 10-3° rs -11 x 10-3°, 
which is in between the above two numbers. 


([Espinoza et al.ll2011ri do not quote a value for 
the decay time scale for this glitch. By expand¬ 
ing and gridding their Figure 5, and reading off 
values of the peak and the first decay point, one 
can obtain an approximate value for r^. This turns 
out to be « 13 ± 3 days, for the decay of both the 
frequency and its derivative (top and bottom pan¬ 
els respectively of that figure). The number de¬ 
rived in this work, = 4.7 ± 0.5 days, is a factor 
of Rs 2.8 smaller, although it is in the right range 


of de cay time scales for Crab pulsar ([Lvne et al 
200011 . If it turns out that the correct value of 


is indeed rs 13 days, then the reason for the factor 
of 2.8 underestimate in this work is most probably 
on account of lack of observations very close to the 
glitch epoch, to which the estimate of Td is very 
sensitive. 


5. Discussion 


In section |4| we concluded that the hard xray 
timing of the Crab pulsar glitch of 15 July 2000, 
using HEXTE/RXTE data, is consistent with the 
results obtained at radio frequencies by the Jo- 
drell Bank observatory. This is, to the best of our 
knowledge, the first time that a glitch of the Crab 
pulsar has been analyzed using xray data, result¬ 
ing in what can be considered to be the closest to 
absolute timing of the Crab pulsar at hard xray 
energies. 


( Rots et al. 200^ found that the main peak of 
the xray pulse profile of the Crab pulsar (the peak 
at phase 0 in Figured]) leads the main peak of the 
radio pulse profil e (just after the ra dio precursor; 
see Figure 6a of ( Kuioer et al. 2003ll ') by 344 ± 40 
^sec. Figure ID shows the result of fitting the 
data of Figure |3| combined with the Jodrell data 
at radio frequencies, obtained from their monthly 
ephemeris. It shows the combined pre and post 
glitch residuals after subtracting the pre and post 
glitch models, respectively. The pre glitch model 
is given in Table 1. The post glitch model is ob¬ 
tained by including additionally in Equation 1, 
permanent step changes in the second and third 
time derivatives of the rotation frequency, Ai/, 
and Az>p, respectively. This solution gave the low¬ 
est rms residual of 3.2 milli periods for the 34 
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post glitch arrival times; the values of the rest 
of the parameters obtained in this ht are consis¬ 
tent with those in Table 2. The mean pre glitch 
separation of the Jodrell data with respect to the 
xray data is 411± 167 fisec. Although this result 
has ss 4.2 times larger error than that quoted by 
( Rots et nil2006l l. it is consistent wit h their result, 
and also with the results of others ( Kuioer et al.l 
2003; Carrillo et al. 2012tl . The post glitch differ¬ 


ences with Jodrell data are difficult to analyze, due 
to the baseline varying in a quasi periodic manner. 
However, after accounting for these variations, the 
last three post glitch Jodrell data certainly sup¬ 
port the above number. 



Fig. 4.— Pre and post glitch xray residuals (this 
work, dots; see text for the corresponding models) 
along with eleven radio residuals from Jodrell Bank 
data (stars), which lie consistently above xray residu¬ 
als, except for the residual of 15 Aug 2000, which has 
a very large error bar (4 ms). 


Several properties of the Crab pulsar can be 
studied as a function of pre and post glitch epochs. 
Three parameters were plotted as a function of 
epoch - (1) The separation of the two peaks in 
the Crab pulsar’s integrated profile, (2) the ratio 
of their peaks, and (3) integrated energy in the 
pulse profile. The first t wo showed no va riation 
worth reporting (see also (iRots et al.ll2006 )). Fig¬ 
ured shows the dead time corrected photon count 
variation of the Crab pulsar as a function of epoch. 

Each point in Figure [5] is obtained by first es¬ 
timating the average off pulse counts in the in- 





^ . 
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Fig. 5.— Normalized and dead time corrected on 
pulse energy of Crab pulsar during the year 2000. 


tegrated profile, then subtracting it from the on 
pulse counts before integrating them, then divid¬ 
ing the result by the average off pulse counts; 
this is done for every DELTAT (16 sec) of data 
for each ObsID. Off pulse counts are obtained 
by integrating in the phase range 0.61 to 0.79 in 
Figure [T] the rest of phase range represents the 
on pulse. Although this procedure will exclude 
any small off pulse em ission from the Crab pulsar 
( Tennant et al.ll20^ . it has the advantage of cor¬ 
recting for dead time, which is not available for 
Barycenter corrected data, and also for data that 
has been filtered in energy range, which is the sit¬ 
uation with our data; see “The XTE Technical 
Appendix” on RXTE websit^Zl- Finally all 48 en¬ 
ergies are normalized by their mean value, which 
lies at the value 1.0 in Figure [3 The assumption 
made here is that dead time correction is, to a 
large order of accuracy, the same for the on and 
off pulse phases of a pulsar, even for bright pulsars 
such as the Crab. Dead time is the duration imme¬ 
diately after the arrival of an xray photon or high 
energy particle, during which the HEXTE detec¬ 
tors are unable to process any more photons. For 
the HEXTE detector there are two sources of dead 
time. One is due to the arrival of an xray photon, 
after which the HEXTE detectors are ’’dead” for 
16 to 30 /Ltsec. The much larger effect is due to 
arrival of high energy particles, after which the 
detectors are dead for 2500 /isec; these are known 
as XULD events. Even for a luminous pulsar such 


^heas arc.gsfc.nasa.gov/docs/xte/appendix_f.html 





























as the Crab and its nebula, the normal photon 
count rate is typically < 400 photons per sec per 
cluster, while the normal XULD event rate is typ¬ 
ically 150 particles per sec per detector (there are 
4 detectors per cluster). Taking the mid value of 
23 fisec, it can be easily seen that the ratio of the 
photon to XULD contribution to the dead time is 
< 0.6%. Since the XULD event rate is the same 
for on and off pulse regions, the above method of 
dead time correction is justified. If required, more 
refined dead time correction can be done by using 
the integrated pulse profile of the pulsar. 

In Figure EJ the standard deviation of the 
spread in values is ~ 16% of the mean value; 
the corresponding spread in the uncorrected fluxes 
is ~ 49% of the mean value. This technique 
of dead time correction for pulsars significantly 
improves the precision of the estimate of inte¬ 
grated xray flux, and may therefore prove use¬ 
ful to find out if the Crab pulsar xray flux varies 
on time scales of decades; for the crab nebula 
this has alr eady been observed a t seve ral xray en¬ 
ergies (see (jWilson-Hodge et al.ll201lh and refer¬ 
ences therein). 

In Figure [5] there is no strong evidence for the 
xray flux of the Crab pulsar to be affected by the 
glitch. The mean pre and post glitch xray fluxes 
are 1.06 and 0.98, respectively, while the standard 
error on these means is 0.03. The mean xray fluxes 
differ by 0.08 ± 0.04, which can not be considered 
a strong result. Moreover, this decrease appears 
to be mainly due to the data well after the glitch 
in Figure El making it less likely to be related to 
the glitch itself. One instrumental feature that can 
simulate such a result is a variation of the average 
off pulse counts as function of epoch in Figure El 
however this does not appear to be the case. Al¬ 
though this can not be considered a strong result, 
it may be interesting to speculate on the possi¬ 
ble causes of glitch related xray flux variations in 
rotation powered pulsars. Such studies have not 
been done so far. Although glitch related xray flux 
enhancemen ts have been reporte d in some Mag- 
netars (see ( Espinoza et al.l[20111 ) and references 
therein), it is generally believed that in Magne- 
tars and AXPs, glitches ar e not always asso ciated 


with xray flux variations ( Dib et al. 20081) . For 


the Crab pulsar glitch under consideration here, 
the post glitch permanent changes in rotation fre¬ 
quency Avp and its derivative Ai>p can contribute 


at most to a fractional increase of ~ 0.01% in the 
post glitch xray flux, while Figure El shows a de¬ 
crease of 8%; so simple energy loss rate of the Crab 
pulsar m ay not be the ex planation. The glitch 
model of ( Ruderman 20091) provides a causal con¬ 
nection between a glitch in a rotation powered pul¬ 
sar and changes in its surface magnetic field. In 
this model, a glitch is caused by excessive stresses 
in the neutron star crust, that are built up by 
super fluid vortices moving outwards due to the 
spin down of the neutron star, dragging with them 
magnetic flux tubes. The crust cracks at the in¬ 
stant of the glitch, and then adjusts itself to a new 
configuration, leading to a corresponding reconfig¬ 
uration of the surface magnetic field, presumably 
both in terms of its field strength as well as in 
terms of its field line structure (curvature of field 
lines, direction of the opening bundle of field lines 
at the surface, presence of higher magnetic multi¬ 
poles, etc). There are several quantitative uncer¬ 
tainties in the predictions from this model. How¬ 
ever it may be worth exploring if a small change in 
surface magnetic field structure of the Crab pulsar 
can lead to a non-linearly large change (decrease 
in the present case, but maybe an increase in other 
glitches) in the post glitch xray flux. 

It is well known that the dispersion measure of 
the Crab pulsar varies d ue to radio propa gation 
within the Crab nebula ( Lvne et al.l Il993 ). and 
that the radio pulses suffer significant and vari¬ 
able refractive and scattering effects within the 
Crab nebula (see ( Backer et al. I [2^ and refer¬ 
ences therein). Such effects can cause errors in 
the arrival times of the radio pulses, with conse¬ 
quent errors on the derived glitch parameters. The 
above interstellar effects are frequency dependent, 
and are absent for pulses at xray energies. The 
consistency between the radio and hard xray tim¬ 
ing results of this work imply that the effect of 
radio pulse propagating through the Crab nebula 
has not been significant during the glitch of July 
2000 . 

Glitches are one of the two timing irregular¬ 
ities observed in rotation powered pulsars, the 
other being timing noise, which manifests obser- 
vationally as random wandering of timing residu¬ 
als. It is currently believed that timing noise is 
due to instabiliti es in the pulsar magn etosphere 
( Lvne et al ] l20inl) : see also (j Arons!l2009l ) and ref¬ 
erences therein for a possible theoretical explana- 
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tion in terms of magnetic reconnection. Now, the 
radio and xray emitting regions in the Crab pul¬ 
sar’s magnetosphere are supposed to be identical 
(the outer gaps), except for the radio precursor, 
which is supposed to arise in the polar gap, while 
for other rotation powered pulsars (eg. Vela) the 
radio and xray e mitting regions are supposed to 
be different; see ( Hardind 2009l l and references 
therein. Therefore simultaneous radio and xray 
studies of timing noise in the Crab pulsar should 
show similar timing noise properties ( rms resid¬ 
uals, time scales of wandering, etc), but probably 
not in Vela like pulsars. 

Glitches, on the other hand, are supposed to 
be due to the steady differential spin down of the 
super flu id core and oute r crust of the neutron 
star; see ( B,uderm^l2009ll and references therein. 
At the instant of the glitch, the outer crust speeds 
up, thereby also speeding up the pulse emitting re¬ 
gions in the magnetosphere, which are firmly an¬ 
chored into the crust by means of the magnetic 
field. Simultaneous timing observations of pulsar 
glitches at radio and xray wavelengths can proba¬ 
bly be used to find out if the emission regions at 
different energies are anchored equally firmly onto 
the surface of the neutron star. For example in 
the case of the Crab pulsar one would expect the 
glitch behavior to be almost identical, while for 
the Vela like pulsars one may or may not notice 
differences. 
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